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The scaling form of the normalized ZFC and FC susceptibility of superparamagnets (SPM's) is presented as a function of the 
normalized temperature y (= kgT /K„(V)), normalized magnetic field h (= H/Hk), and the width cr of the log-normal distribution 
of the volumes of nanoparticles, based on the superparamagnetic blocking model with no interaction between the nanoparticles. 
Here (V) is the average volume, K u is the anisotropy energy, and H K is the anisotropy field. Main features of the experimental 
results reported in many SPM's can be well explained in terms of the present model. The normalized FC susceptibility increases 
monotonically increases as the normalized temperature y decreases. The normalized ZFC susceptibility exhibits a peak at the 
normalized blocking temperature yi, (= kBTb/K u (V)), forming the yj vs h diagram. For large cr (cr > 0.4), y& starts to increase 
with increasing h, showing a peak at h — hb, and decreases with further increasing h. The maximum of yb at h - hb is due to the 
nonlinearity of the Langevin function. For small cr, yj monotonically decreases with increasing h. The derivative of the normalized 
FC magnetization with respect to h shows a peak at h — for small y. This is closely related to the pinched form of Mfc vs H curve 
around H — observed in SPM's. 
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It is well known that the magnetic susceptibility of super- 
paramagnets (SPM's) and spin glasses (SG's) strongly depends 
on the way how the cooling of the system is done before or 
during the measurement [1]. Typically there are two types of 
susceptibility, zero-field-cooled (ZFC) susceptibility and field- 
cooled (FC) susceptibility. In the case of the ZFC susceptibility 
CtzfcX the system is cooled in the absence of an external mag- 
netic field (H) from high temperatures well above a freezing 
temperature Tf (a blocking temperature T/, for the SPM's and 
a spin freezing temperature T g for the SG's) to the lowest tem- 
perature well below Tf (the ZFC cooling protocol). After H 
is applied at the lowest temperature, the ZFC susceptibility is 
measured with increasing temperature (T). In contrast, the FC 
susceptibility (xfc) is measured in the presence of H with de- 
creasing T from high temperatures well above Tf to the lowest 
temperature well below Tf (the FC cooling protocol). 

It has been experimentally confirmed that the T dependence 
of the susceptibility for the SPM's is essentially different from 
that for the SG's E Ql HI El [6] El El 13 [JOl |TT] [121 . For SPM's 0, 
the ZFC susceptibility exhibits a peak at the blocking temper- 
ature T/,, while the FC susceptibility monotonically increases 
with decreasing T, showing no anomaly at T\,. The difference 
between the FC and ZFC susceptibility appears well above Tb- 
For SG's [3 1, the ZFC susceptibility also exhibits a peak at the 
spin freezing temperature T g , while the FC susceptibility in- 
creases with decreasing T above T„ and starts to decrease or to 
saturate below T g , indicating the evidence of glass transition. 
The difference between the FC and ZFC susceptibility appears 



just at T g . 

The purpose of the present paper is to explain the peculiar 
features on the T and H dependence of the ZFC susceptibil- 
ity, which were reported in many SPM's [3] S |5j |7] [8] 
13 [10] El] HH. The first feature is as follows. The ZFC sus- 
ceptibility shows a peak at the blocking temperature T/,. The 
blocking temperature Tb changes with H in two different ways, 
depending on the nature of SPM's. (i) Tb first increases with 
increasing H, reaching a peak and then decreases with further 
increasing H (unusual type), (ii) Tb monotonically decreases 
with increasing H (conventional type). The feature (i) is ob- 
served for SPM's such as Fe304 nanoparticles (Luo et al. |2]), 
natural horse-spleen (Friedman et al. [7]), y-Fe203 nanoparti- 
cles with diameter 7 nm and ferritin (Sappey et al. El), and a 
diluted magnetic fluid composed of FePt nanoparticles (Zheng 
et al. Ifl2l ). These two different features are observed even in 
the same SPM, depending on the conditions of samples (such 
as dilute and concentrated samples). Fory-Fe203 nanoparticles 
in dispersed polymer show that the feature (i) is observed for 
dilute samples and that the feature (ii) is observed for concen- 
trated samples (Kachkachi et al. IfTTI ). The second feature is the 
H dependence of the FC magnetization Mfc near H — 0. The 
derivative of Mfc with respect to H has a sharp peak at H = 
for ferritin (Tejada et al. |9|) and for CoFe204 nanoparticles 
embedded in potassium silicate (Zhang et al. 0). For natu- 
ral horse-spleen the magnetic hysteresis shows an anomalous 
pinched hysteresis loop near zero field (Friedman et al. Q). 

Here we present our numerical analysis of FC and ZFC sus- 
ceptibility, based on the superparamagnetic blocking model 
with no interaction between the nanoparticles. The theory used 
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here is the same as that used by Bitto et al. J3j, and is essentially 
the same as the Neel-Brown model |[T3l[T4l . We show that the 
main features of the experimental results of the ZFC and FC 
susceptibility for SPM's can be well explained in terms of the 
present model. The peculiar features of the susceptibility is due 
to the Langevin function of the magnetization. Our numerical 
results will be compared with those predicted by Sappey et al. 
1 8 1 . We show that the blocking temperature Tb is extremely sen- 
sitive to the width cr in the log-norm distribution for volumes of 
nanoparticles. Resonant tunneling is one of the models to ex- 
plain the shift of Tb with H l6l [T5ll . For example, 7j is 65 K 
for y-Fe2C>3 (Sappy et al. |8|) and is too high for the quantum 
tunneling by which the magnetic moment of the particles can 
flip. So this possibility may be ruled out. 

2. Scaling form of Xzfc and \fc 

2.1. Magnetization in the SPM state 

We consider a single-domain particle (nanoparticle) with vol- 
ume V and the anisotropy having uniaxial symmetry. We as- 
sume that there is no interaction between these domains. The 
simplest form of the anisotropy energy is given by K U V sin 2 0, 
where is the angle between the magnetization direction and 
the easy axis, and K u is the anisotropy energy per unit volume. 
The two states [0 = (partallel) and <p = n (antiparallel)] as 
the ground state, are energetically degenerate and separated by 
a energy-barrier height AE CI {V) = K U V. At temperatures well 
above Tb, the magnetization of the single-domain fluctuate due 
to the thermal activation energy (kgT) in a paramagnetic way 



(SPM state). The definition of Tb will be given in Sec. 2.3 The 
magnetization (the magnetic moment per unit volume) in the 
SPM state is described by 



with 
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M spm = Mj ( cos 0) ; 

JV^ cos 6dQ 
Je'^dQ, 

n r M S VH cos B 

J e v cos 0(2n sin 6dG) 
o 

^ M S VH cos B 

J e v (2n sin 8d6) 



k B T 



(1) 



(2) 



where // (= M S V) is the magnetic moment of the single do- 
main with the volume V, M s is the saturation magnetization 
(the magnetic moment per unit volume), (0 < 6 < n) is 
the angle between // and H (the z axis), Q is the solid angle, 
dQ. = 2n sin 9d9, ks is the Boltzmann constant, U is the Zee- 
man energy defined by 



U = -fi-H = -M s VHcost 
and L(£) is the Langevin function of (; 

L(0 = coth(<D - I 



(3) 
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2.2. Magnetization in the blocked state 

At temperatures well below Tb, this themal activation energy 
is too small to cause the barrier hopping between the two states 
(0 = and n). Then the direction of the magnetization in the 
single domain is frozen to one of the two states (0 = and 
4> = n). 



2.2.1. H || easy axis 

Here we evaluate the ZFC susceptibility below Tb in the sys- 
tem where the single domains are randomly oriented ifTBI IT71 
Q~8l[T9). First, we consider the case when an external magnetic 
field H is applied along the easy axis. The system is in ther- 
mal equilibrium. The total energy F\\ consists of the anisotropy 
energy and the Zeeman energy; 



F» = V(K U sin 2 



HM, cos i 



K u V(sin z 



and the derivative is given by 
dF« 



2K u V{cos <p + h) sm ' 



2/icos0), (5) 



(6) 



where <p is the angle between the easy axis and the magnetiza- 
tion direction. The normalized field h is defined by h - H/Hk, 
where Hk is the anisotropy field and is defined by 



H K = 2KJM S 



(7) 



The free energy F\\ has a minimum at = 0, and n, and a 
maximum at cos = —h. The minimum of the free energy F\\ is 
given by 

F™° = K u V(-2h), (8) 

at = and n, while the maximum of the free energy is given 
by 



F,7 X = K u V(h 2 + 1). 



(9) 



at cos = —h. The energy-height barrier between two states 
(0 = and ji) is obtained as 



AF,| =F™ x -K™ = ^V(l+/z) z 



(10) 



This means that the energy-barrier height increases with in- 
creasing h. Then the magnetization along the easy axis (0 = 
state) is simply given by 



My = M s . 



(ID 



The appearance of magnetization along the easy axis below Tb 
is indicative of the freezing of the direction of the magneization 
along the easy axis. This is one of the important features for 
the blocked state. The susceptibility along the easy direction is 
equal to zero; 

X\\ = 0. (12) 
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2.2.2. H ± easy axis 

Next we consider the case when the magnetic field is applied 
to the hard axis (perpendicular to the easy axis). The total free 
energy F ± is given by 



F ± = V(K U sin 2 4> - HM S sin 0) = K u V(sm 2 , 



and the derivative is given by 
dF ± 



• - 2h sin i 



dcf> 



2K u Vcos <p(s,in<p- h), 



(13) 



(14) 



where <p is the angle between the easy axis and the magnetiza- 
tion direction, and h is smaller than 1 . The free energy F± has 
a minimum at sin0 = h. Then the magnetization along the hard 
axis (the field direction) of the single domain is obtained as 



M s 

M ± = M s sin <f> - M s h = — -H, 
Hk 



(15) 



while the susceptibility along the hard axis (the field direction) 
is 

_ M x _ M s _ M] 
X± ~ ~H ' " ~H~ K ~ 2K U ' 
The minimum value of F x is given by 



(16) 



Hcos e 




7L ± H sin 8 
H sin 9 * 

Figure 1 : (Color online) The schematic diagram of the easy axis (the z axis) of 
the single domain and the external magnetic field H in the z-x plane. The angle 
between H and the z axis is 8. The z and x components of H are H cos 8 and 
H sin 9, respectively. Correspondingly, the z and x components of the magneti- 
zation induced by the field H are^yff cos 6 and x±H sin 8. From the symmetry, 
only the component of the resultant susceptibility along the direction of H is 
not equal to zero, and is given by x\\ cos 2 8 + x± sin 2 8. 



F'l m = -K u Vh 2 , (17) 

at sin 4> - h, while the maximum value of F ± is given by 

FT* = K u V(l -2/0, (18) 

at <p — n/2. Then the energy-barrier height between two states 
denoted by sin 9 - h, is given by the energy difference AF ± 
defined by 



AFj_ = F™ ax - Ff n = AE a (V, h) = K u V(l - h) 1 



(19) 



indicating that the energy-barrier height AE a (V,h) decreases 
with the increase of h or with the decrease of V. 

In a system where the direction of the easy axes for single 
domains is randomly oriented, the average susceptibility (so- 
called powder susceptibility) can be evaluated as 



Xa 



^ll<cos 2 6»>+^ ± <sin 2 6») 
JOtfl cos 2 6 +x±. sin 2 6)dQ. 



JdCl 



1 2 

3*1 + ^ 

M] 



(20) 



using Eqs.([T2J and (JT6) , where dQ. = 2n sin 8d8, 6 (0 < 6 < n) 
is the polar angle between the easy axis (the z axis) and the 
magnetic-field direction, see Fig. [T] and the caption for the de- 
tail of calculation. When the average susceptibility is expressed 



by Eq.(20i, we refer to the system as being in the blocked 



state. This susceptibility is the ZFC susceptibility of the single- 
domain with volume V below T\,. 



2.3. Blocking temperature Tb 

The relaxation time of the magnetization between two states 
(0 = and n states) is given by thermal activation (Arrehenius 
law). In the Neel-Brown relaxation process lfT3lfT4l . the relax- 
ation time of the magnetization between the two states is given 

by, 

AF.JV) 

(21) 



r exp[ ; - ], 



k B T 

where to is a microscopic limiting relaxation time (usually to ~ 
10~ 9 sec). The measurement time r m is typically of the order of 
10 2 sec for the DC magnetization measurement. In Eq. (21 1 
with t — T m and T = Tt,(V), the blocking temperature Tb(V) is 
derived as 

AE a (V) 



T h (V) = 



(22) 



k B ln(T m /T ) 

In the presence of H, it is assumed that the energy-barrier 
AE a (V, h) is described by 



AE a (V,h) = K u V(l-hT, 



(23) 



where a — 2 in the present case. Note that a — 1.5 when 
the orientational order is taken into account [8]. The blocking 
temperature Tb(V, h) is given by 



T b {V,h) 



K U V 



k B ln{T m /T ) 



(1-/0" 



(24) 



2.4. Form of magnetization in the superparamagnetic state and 
blocked state 

Before discussing the general case, first we consider a simple 
case where all nanoparticles have the same volume. Above the 
blocking temperature !/,( V, h), the system is in the SPM state. 
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Figure 2: (Color online) Plot of f(x,cr) as a function of x at various cr (cr = 
0.1 - 0.8). f(x,<r) is a log-normal distribution function and is denned by Eq. 



The magnetization M 1 '*"" (per unit volume) for the SPM state is 
given by 



M s z p F m c (V) = M^ m (V) = eM s L( 



FC 



M S VH 



(25) 



where e is the volume fraction occupied by ferromagnetic 
nanoparticles. 

Below Ti,(V,h), the system is in a blocked state. The mag- 
netization M bl (per unit volume) of the blocked state is given 
by 



M°'(V) 



3K„ 



(26) 



and 



M F ' C (V) = eM s L( M /*,J . 



k B T b (V,h) 

2.5. Scaling form of normalized ZFC and FC susceptibility 



(27) 



Using Eqs. (25i-(27i, we derive the general expression for 
the magnetization. We consider a system of ferromagnetic 
nanoparticles having a wide distribution of volume sizes and 
let (V) be their average volume. We also define a characteristic 
volume V m (T, h) given by 



V n (T,h) 



knT 



K u (l - Kf 
which is derived from the condition 

AE a (V m ,h) 



ln(T m /r ), 



k B T 



ln(T m /T ). 



(28) 



(29) 



For simplicity we introduce a volume ratio x (= V/(V}). The 
characteristic volume ratio x m is then defined by 



V m (T,h) 



knT 



(V) 

y 



(1 - h) a 



K u (V)(l-h)<* 
ln(T m /T ), 



MTm/To) 



(30) 



where y is the reduced temperature and is defined by y — 
k B T/(K u {V)). For x < x,„, the system is in a SPM state, and for 
x > x m , the system is in a blocked state. For the SPM state, the 



ZFC and FC magnetizations shown in Eq. ( 25 1 can be rewritten 

as 

2hx 

' : (x, y, h) = M s F P ™(x, y, h) = eM s L{ — ), (31) 



M spm ( 



while for the blocked state Eqs. ( 26 1 and ( 27 1 are rewritten as 

2/7 



M° FC (x,y,h) 



3 ' 



and 



M F ' c (x,y,h) = eM s L( 



2h 

(i - hy 



ln(T m /r )), 



(32) 



(33) 



respectively. 

We assume that volume distribution of the nanoparticles, 
f(x, cr), is expressed by the log-normal function [3 1 



f(x,a) = ——— exp[ — —-} 
S1ko-x 2cr z 



(34) 



where 



lA. 1 

/ 



f(x, cr)dx = 1 , and x a 



s 



xf(x,cr)dx = exp(cr /2). 



Here x av is the average of x and cr is the width of the distri- 
bution. Figure [2] shows a plot of f(x, cr) as a function of x at 
various cr. The log-normal distribution function f(x, cr) has a 
maximum at x — x max = exp(-cr 2 ). The value of x max de- 
creases with increasing cr, while the value of x m , increases with 
increasing cr. In the limit of cr — > 0, f(x, cr) becomes a Dirac 
delta function 5{x - 1) which has a sharp peak at x — 1. Using 



Eqs. (31 i-(34i, the scaling forms of the normalized ZFC and 
FC susceptibility are given by 

XzFc(y,h,cr) 1 r°° r ,.2/ix 

[L( )U-i(x m -x) 

Jo 



Xo 



h Jo y 

2h 

+ —U- l (x-x m )]f(x,cr)dx, (35) 



and 



XFc(y, h, cr) 1 r°° 2hx 

- — = 7 [L( )U- X {x m -x) 

Xa « Jo y 



,, o; ln(Tm/To) 

L(2n— ^-)U-i(X — x m )lf(x, cr)dx, 



(1 - h) a 



(36) 



where U-\ is a step function [U-\(x) = 1 for x > and for x < 
0]- Xo i- /(2K U )] is a constant susceptibility and x m given 



by Eq. (30 1 is a function of y, h, and t,„/to. The normalized 
ZFC and FC susceptibilities depend only on y, h, T m /ro, and 
cr. They are independent of the values of e, M s , K u , and (V). 
Hereafter, we assume ln(r m /ro) = ln(10 2 /T0~ 9 ) = 25.328 for 
simplicity. In this case, Xzfc/Xo ^ndxFc/Xo aie described by 
scaling functions of y, h, and cr. 



3. Numerical calculation 
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Figure 3: (Color online) Plot of Xzfc Ixo . Xfc/Xo, and S = (xfc ~Xzfc)Ixo 
as a function of reduced temperature y. h = 0.1. The width tr is varied from 
0.1 to 0.6 as a parameter, a = 2.0. ln(r m /r ) = 25.328. xo = eM 2 J(2K„), 
y = k B T/(K u (V)), and h = H/H K . 
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Figure 4: (Color online) Plot of Xzfc/Xo, Xfc /Xo , and & = (Xfc ~ Xzfc)/X0 
as a function of reduced temperature y. a = 0.6. The reduced field h is changed 
as a parameter: h = 0.005 - 0.5. a = 2.0 



Figures |5Ja)-(f) show typical plots of Xzfc/Xo, Xfc/Xo and 
the difference 6 defined by 5 — (xzfc ~Xfc)/Xo as a function of 
the normalized temperature (y = kgT / K U (V)) , where a = 2.0 
and ln(r m /To) = 25.328. The reduced field h (= H/Hk) is kept 
constant at h = 0.1, and the width of the log-normal volume 
distribution function (cr) is varied as a parameter (cr = 0.1 to 
0.6). The normalized ZFC susceptibility Xzfc/Xo f° r cr = 0.1 
exhibits a peak at the normalized blocking temperature y^,. The 
difference 5 starts to appear at the onset normalized temperature 
of irreversibility y !rr (just above yb) and increases with decreas- 
ing y. The irreversible effect of susceptibility occurs below y lrr . 
The normalized FC susceptibility Xfc/Xo f° r °" = 0.1 is nearly 
independent of y below yb- These features are common to those 
observed in real SG systems J3J. Here, the system behaves like 
a superspin glass (SSG), where the large magnetic moment of 
the nanoparticle plays the same role as an atomic spin in a SG 
system l20l . As cr increases, the peak of Xzfc/Xo becomes 
broad and shifts to the large y-side. The normalized FC suscep- 
tibility Xfc/Xo for large cr (typically cr = 0.6) tends to increase 
with decreasing y in the small-y region (at low temperatures). 

Similar numerical calculations are carried out for the differ- 
ent a (= 1.5-2.0). We find that the essential points of numerical 
results are not so different for different a. So here we show only 
the result for a = 2.0. Note that a — 1.5 is used in the paper of 
Sappey et al. (5). 

Figures|4](a)-(f) show typical plots of Xzfc/Xo, Xfc /Xo and 6 
as a function of y for cr = 0.6, where h is changed as a parameter 
(h = 0.005 - 0.5). The normalized susceptibility Xfc/Xo at 
each h increases with decreasing y below y^. The normalized 
susceptibility Xzfc/Xo exhibits a broad peak at yb, which shifts 
to the low-y side with increasing h. 

Figures [5ja) and (b) show a plot of yb as a function of h 
(yb vs h diagram) at various width cr. For small cr (typically 
cr < 0.11), yb monotonically decreases with increasing h. At a 
larger cr (typically 0.6 < cr < 1 .2), yb decreases with increase in 
h for very low h, and shows a local minimum around h - 0.05. 
With further increasing h, yb increases with increasing h, show- 
ing a local maximum around h - 0.10 - 0.15, and decreases 
with further increasing h. The local maximum position shifts to 
high-/; side with increasing cr. The origin of the local maximum 
in the yb vs h curves arises mainly due to the nonlinearity of the 
Langevin function in Eq. (33 1, as will be discussed in Sec. [4] 
It should be noted that this diagram of yb vs h may be a fea- 
ture common to any SPM's, since it depends only on the scaled 
variables h and y. 

As shown in Fig.|5Ja), the curvature of yb vs h diagram drasti- 
cally changes with varying cr. We assume that the h dependence 
of yb is expressed by a power law form 



h = h (l 



yb 
yo' 



f. 



(37) 



for 0.1 < h < 0.6, where /3 is an exponent and ho and yo are 
constants. The least-squares fit of the numerical result of h vs 
yb for each cr to Eq. (37 1 yields the parameters yo, and ho. 



In Fig. [6] we show the value of (3 thus obtained, as a function 
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Figure 5: (Color online) (a) yy vs h diagram at various cr (= 0.05 - 1.2). yj, is 
the normalized blocking temperature, and h is the normalized field, a = 2.0. 
ln(r„,/To) = 25.328. cr is the width of the log-norm distribution, (b) Detail 
of the yi, vs h diagram at a = 0.6, 0.7, 0.8, 1.0, and 1.2. The dotted lines are 
denoted by yt = 2h/(o with £q = 1, 2, • • • , 8. 
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Figure 6: (Color online) Exponent ji vs cr. The exponent p is derived from the 
least-squares fit of the yj, vs h curves for 0.1 < h < 0.6 to Eq. J37|. a = 2.0. 



Figure 7: (Color online) (a) y b (cr, h) vs cr at h = 0.001, 0.1, 0.2, 0.3, and 0.4. a 
= 2.0. (b) Plot of = yi,(<r,h)/yi,((T — > 0, h) vs cr at various h. yb(cr,h) is the 
normalized blocking temperature at fixed h and tr. 



of cr. We find that the exponent j3 is strongly dependent on the 
width cr: ft increases with decrease in cr. The exponent f3 is 
nearly equal to 0.5 at cr — 1.2. It increases with decreasing cr 
and is nearly equal to 1.5 at cr = 0.01. Note that [3 — 3/2 is an 
exponent predicted for the de Almeida-Thouless (AT) line [21 J 
for the H-T diagram in Ising SG systems. 

Figure |7|a) shows the plot of yb(cr, h) as a function of cr for a 
fixed h(- 0.001-0.4). The value of yt, increases with increasing 
cr at each h. Figure |7|b) shows the ratio fi/, defined as = 
yh(cr, li)/yb(cr — > 0, h) as a function of cr at various h, where 
yb(c — » 0, h) is obtained from the extrapolation of yb(cr, li) in 
the limit of cr — > 0. We find that the ratio fit, is proportional 
to cr at small cr and tends to saturate at large cr. Our result is 
rather different from the prediction by Sappey et al. [8 1 above 
cr = 0.8. Their value of tends to diverge with increasing cr 
above cr = 0.8, while that of ours tends to saturate. The reason 
for such a difference in ptb at large cr is that the nonlinearity of 
the Langevin function with y and h is not taken into account in 
the paper by Sappey et al. O. 

Figure [8ja) shows the h dependence of the derivative 
d(hxFc/Xo)/dh at various y, where cr = 0.6, a = 2.0, and 
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Figure 8: (Color online) /i-dependence of (a) d(hxFc/Xa)/dh and (b) Figure 9: (Color online) (a) The peak height ij p and (b) the peak field h p as a 
d(hXzFc/Xo)/dli at various y (= 0.01 - 0.5). cr = 0.6. a = 2.0. yt is nearly function of y at various cr, where AQixzfc IXo)l&h vs h curve exhibits a peak 
constant (as 0.065) for 0.02 < ft < 0.12 (seeFig.|5j. height^,, at a peak field h p for each y. a = 2.0. 
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QiXfc/Xo) is proportional to M FC . The value of yj is nealy 
constant (» 0.065) for 0.02 < /i < 0.12 (see Fig.|5j>. For y < y fc , 
it shows a sharp peak around h — 0. For y > y^, d(hxFc/Xo)/dh 
becomes fiat around h = 0. The pinched hysteresis form for 
M FC vs H near // = observed for SPM's EdlSl implies that 
the derivative of M FC with respect to H exhibits a sharp peak. 
Thus these experimental results can be explained by the present 
model. Figure [HJb) shows the h dependence of the derivative 
d(hXzFc/Xo)/dh at various y, where cr = 0.6 and (hxzFc/Xo) 
is proportional to Mzfc- It shows a broad peak at h — h p (= 
0.0159) for y = 0.03. This peak field h p shifts to the lower- 
h side with increasing y, while the peak height rj p increases 
with increasing y, showing a local maximum at h p = 0.0089 
for y w 0.065. We note that these values of y and h p for the 
local maximum are located on the yj, vs h diagram (h - 0) for 
cr = 0.6 (see Fig. |5|b)). Figures |9ja) and (b) show the peak 
height rj p and the peak field h p as a function of y at various 
cr for d(hxzFC /Xo)/dh vs ^> respectively. The peak height j] p 
shows a local maximum around y = 0.05 - 0.06. The peak field 
rj p decreases with increasing y. These features are independent 

of CT. 

In summary, we obtain the following important features from 
the above numerical calculations, (i) The normalized ZFC sus- 
ceptibility XzfcIxo vs y nas a broad peak at y = y&, forming 
the yt, vs h diagram. The value of yt, exhibits a local minimum 
and a local maximum for large cr (cr > 0.4). (ii) The deriva- 
tive &(hxFclXo)l&h vs h shows a very sharp peak around h = 
at low y. (iii) The derivative d(hxzFC /Xo)/dh vs h shows a 
peak at the boundary of the y& vs h diagram (h 0). (iv) The 
monotonic increase of Xfc/Xo with decreasing y is seen below 
yt, for relatively large cr. With decreasing cr, Xfc/Xo becomes 
flat below yt, like SG's. 



4. Discussion 

Both the non-monotonic H dependence for the ZFC-peak 
temperature Tt (2] |8] [TT] [12) and the sharp peak of the 
derivative dM/dH have been observed in many SPM's 00|9). 
Several theories have been presented for the explanation of 
the maximum of 7$ at low H (8] [TT] [TBI including the theo- 
ries based on the Neel-Brown model |[T3l [T4l and the resonant 
spin tunneling theory ||6] [T3) . The peak of Tt, at small H in 
Mni 2 Oi2(CH3COO)i 6 (H 2 0)4 (denoted as Mn 12 ) @ and ferritin 
|9| may be explained in terms of thermally assisted, field-tuned 
resonant tunneling between particles. However, this model may 
not be valid for the explanation of similar behaviors in SPM's 
with Tt, which is too high for the quantum effect to appear. 

In Sec. [3] we show that these features can be well explained 
in terms of the scaling form of Xzfc/Xo w&XfcIxo- The local 
minimum and local maximum of yt, vs h tends to disappear as 
cr becomes smaller. The non-monotonic behavior of y& vs h is 
mainly due to the non-linearity of the Langevin function. The 
curvature of yj vs h is strongly dependent on cr. It is interesting 
to discuss where the nonlinearity of the Langevin function is 
significant in the yt, vs h diagram. Note thai xzfc I Xo is given 



by Eq. ( 35 i using the Langevin function which is a function of 

f = (M s V)H/(k B T) = 2hx/y, (38) 

where x = V/<V>, y = k B T/(K u (V)), h = H/H k , and H k = 
2K U /M S . M S V is the magnetic moment of the particle of volume 
V and M s is its magnetization per unit volume. The nonlinearity 
is considered to appear when f > 1 . For convenience, we as- 
sume that x = 1, which means V = (V). Then f is approximated 
as £ = £o = 2h/y. In Fig. [5ja), we make a plot of the straight 
lines denoted by the relation £b = 2/j/y, where £o = 1-8. It 
is clear that the local maximum and local minimum for each cr 
are located on the lines with fo = 1 an d fb = 3-4, respectively. 
This implies that the nonlinearity of the Langevin function is 
crucial to the occurrence of the local maximum and local mini- 
mum. In the numerical calculation, Sappey et al. [8 1 have used 
the approximation of the Langevin function for f «: 1, which 
may lead to results inconsistent with our numerical results. We 
note that yt, with cr > 0.4 exhibits a local minimum around 
h — as shown in Fig. |5|a). This local minimum is located 
around the line given by £o = 1 in the y& vs h curve with each 
cr. Experimentally Luis et al. [10] have reported that in natural 
horse-spleen ferritin the curve of Tt, vs H exhibits a local mini- 
mum at H — 0.5 kOe and a maximum at H — 3 kOe. This result 
is qualitatively consistent with the results of our numerical cal- 
culation. 

5. Conclusion 

We show that remarkable features in the T and H dependence 
of the ZFC and FC susceptibility for SPM's (nanopraticles) can 
be well explained by the present model. These behaviors arises 
mainly from the nonlinearity of Langevin function, but not due 
to the quantum tunneling effect. 
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